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In the next generation of long-baseline neutrino oscillation experiments, aiming to determine the 
charge-parity violating phase Scp in the appearance channel, fine-grained time-projection chambers 
are expected to play an important role. In this Letter, we analyze an influence of realistic detector 
capabilities on the Scp sensitivity for a setup similar to that of the Deep Underground Neutrino 
Experiment. We find that the effect of the missing energy, carried out by undetected particles, is 
sizable. Although the reconstructed neutrino energy can be corrected for the missing energy, the 
accuracy of such procedure has to exceed 20%, to avoid a sizable bias in the extracted Scp value. 

PACS numbers: 14.60.Pq, 14.60.Lm, 13.15.+g, 25.30.Pt 


The matter-antimatter asymmetry in the Universe is 
an outstanding problem of modern physics. It is expected 
that equal amounts of matter and antimatter were pro¬ 
duced in the Big Bang, yet we observe a baryon asymme¬ 
try of the order HU 10 . This requires a dynamic mecha¬ 
nism for baryogenesis, a prerequisite for which is violation 
of charge-parity (CP) symmetry p]. While the contribu¬ 
tion of the quark sector is too small by several orders of 
magnitude [23, leptogenesis offers a viable alternative to 
generate the asymmetry (3j. 

Under the assumption of three-neutrino mixing and 
Majorana masses, possible sources of CP violation in the 
lepton sector are the Dirac CP phase 5c p, testable in 
neutrino oscillation measurements, and the Majorana CP 
phases, entering lepton-number violating processes only. 
Owing to the large value of the $13 neutrino-mixing an¬ 
gle [IHB], the 5cp phase has the potential to give an im¬ 
portant, or even dominant, contribution to the matter- 
antimatter asymmetry of the Universe [7]. 

To observe CP violation in neutrino oscillations, an ap¬ 
pearance experiment is necessary [8j. In this Letter, we 
consider the — > i/ e and —> v e transitions, for which 
the oscillation probabilities in vacuum can be approxi¬ 
mated by jlj 

P^ e ~ S23 sin 2 2013 sin 2 A 31 + C23 sin 2 20 12 sin 2 A 2 i 
+ J cos (fpScp A31 ) A 2 i sin A 31 , (1) 

with Sij = sin 0^, dj = cos 0 y , A= (m 2 —m 2 ) L/(4:E U ), 
and J = C13 sin 20 i 3 sin 20 i2 sin 20 23 . Here, E u is the neu¬ 
trino energy and L denotes the distance from the neu¬ 
trino source to the detector. Finally, the upper (lower) 
sign of 5cp refers to the neutrino (antineutrino) chan¬ 
nel. Note that in our calculations, the exact formulas 
including matter effects are used instead of Eq. ([I]). 

Unless 5cp = 0 or tt , the CP symmetry is violated and 
the oscillation probabilities (|T]) are different for neutrinos 
and antineutrinos. For a maximally CP-violating value 
of 5c p, a combination of ongoing oscillation experiments 


such as T2K [TO] and NOvA [TT] will probe CP viola¬ 
tion at the ~ 2cr confidence level. Interestingly, recent 
results from the T2K experiment combined with reactor 
measurements for 0i3 show some preference for maximal 
CP violation, 5cp = —90° [12, [T3]. Nevertheless, the 
accuracy of T2K and NOvA in determining the value of 
5c p will be between 30 and 70° at the lcr confidence 
level BUS]- 

An accurate measurement of 5cp has important conse¬ 
quences for model building. In many flavor models, par¬ 
ticular relations—called sum rules—take place between 
the mixing angles and 5cp HO BE]. Furthermore, in 
models of leptogenesis with sizable flavor effects, the val¬ 
ues of the oscillation parameters have a strong impact on 
the matter-antimatter asymmetry of the Universe mm- 
Therefore a precise determination of 5c p, together with 
the mixing angles and light neutrino masses, can be a 
powerful tool to discriminate among different models and 
help to shed some light on some of the open questions of 
the Standard Model f20l [2Tj . 

A precise determination of the value of 5cp in future 
experimental programs [221425] will require to achieve an 
unprecedented accuracy, keeping systematic uncertain¬ 
ties under control at the percent level. Since neutrino 
beams are rather broad in energy, for any given event 
observed at the detector the neutrino energy needs to 
be reconstructed from the measured kinematics of the 
measured particles in the final state. Because the CP- 
violating phase enters the oscillation probabilities ([I]) 
with a non-trivial dependence on the neutrino energy, 
a bias in the energy reconstruction translates into a bias 
in the determined 5cp value, as we discuss in this Letter. 

Provided the energies of all the particles produced in 
the event are measured, the neutrino energy in charged- 
current (CC) processes can be simply determined using 
the calorimetric method, 

E°*=E t + Y i T” + e n + Y i E j , ( 2 ) 

* 3 
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summing the charged lepton energy Eg, the kinetic ener¬ 
gies of the knocked-out nucleons JA, the corresponding 
separation energy e n , and the total energy of any other 
particle produced Ej. 

Unlike the kinematic energy reconstruction, based on 
the charged lepton kinematics only—typically used for 
quasielastic events |2fif[2t)] -the calorimetric method ([2]) 
is applicable to any final state. This feature is of great 
importance for oscillation studies performed in the high- 
energy regime, where quasielastic processes give only a 
small contribution to the inclusive CC cross section m- 
One of the main advantages of a fine-grained time pro¬ 
jection chamber (TPC) with respect to those detectors 
which only can track the leading charged lepton is its 
ability to distinguish u e CC events from neutral-current 
(NC) backgrounds, even for CC events which are ei¬ 
ther resonant or deeply inelastic. In these types of 
events, multiple tracks can occur and the identification of 
the leading electron track requires fine-grained detection 
technique. Also, this provides the ability to measure the 
energy deposited in the hadron shower. 

Moreover, for an appearance experiment with neutrino 
energies in the GeV range, a fine-grained TPC would in 
principle outperform a much larger water-Cherenkov de¬ 
tector thanks to the ability to reconstruct also events 
which are not quasielastic. However, an accurate recon¬ 
struction of hadrons is a formidable experimental task, 
especially in the case of multi-track events. One needs to 
keep in mind that neutrons typically escape detection and 
any undetected pion leads to an energy underestimation 
by at least its mass, ^135 MeV. The limited detection ef¬ 
ficiencies for the different hadrons produced in the events 
will also contribute to the missing energy budget. Note¬ 
worthy, the number and energy distribution of hadrons 
in general (and neutrons in particular) is very different 
for neutrino and antineutrino CC events. 

In principle, part of these effects may be alleviated 
by using near-detector data and by an accurate determi¬ 
nation of the detector response to different test beams. 
Nevertheless, the near and far detectors will most gener¬ 
ally not be identical in design or in performance, which 
leads to notable uncertainties when determining the de¬ 
tector capabilities. In particular, different dimensions 
of the near and far detectors may result in a signifi¬ 
cantly different containment of neutral pions and neu¬ 
trons. Therefore, the effect of neutral secondaries, like 
neutrons, has to be corrected for by using the detector 
Monte Carlo, which ultimately relies on an event gener¬ 
ator. If the physics model in the event generator pro¬ 
vided an accurate description of the underlying physics 
this would not present a major problem. However, cur¬ 
rently, event generators may not be able to provide suf¬ 
ficiently accurate predictions for the multiplicities and 
energy distributions of neutral secondaries and absorbed 
pions. While in the coming years, more sophisticated 
theoretical models can be expected to be implemented in 


generators, their accuracy would have to be validated by 
data and cannot be relied on a priori. The spirit of this 
Letter is, first, to demonstrate that there is a problem 
arising from missing energy and, second, to explore the 
level of accuracy required in estimating this missing en¬ 
ergy to avoid a deleterious effect on the measurement of 
the CP phase. However, we are not concerned with how 
this can be practically achieved. The detailed simulation 
and study of the effect of the near detector on the de¬ 
termination of missing energy in neutrino interactions is 
well beyond the scope of this work, and will eventually 
have to be performed by the experimental collaborations. 
In this work, we will instead demonstrate and quantify, 
from a phenomenological point of view, how an under¬ 
estimation of the missing energy in neutrino events may 
affect the extraction of the value of 8c p- 

To analyze the effect of realistic detection capabili¬ 
ties on the energy reconstruction in a fine-grained TPC, 
we take into account energy resolutions, efficiencies, and 
thresholds for particle detection. For all hadrons, we set 
them to optimistic values, detailed in Ref. pTTj . In our 
considerations, finite detector resolutions smear particle 
energies according to the normal distributions centered 
at true energy values. 

We make the assumption that all neutrons escape de¬ 
tection. It should be stressed that this assumption is 
rather conservative. For instance, in Ref. [SJj it is esti¬ 
mated that only a 10% of the neutrons with energies be¬ 
low a GeV will escape detection in a liquid argon TPC. 
However, since neutrons travel some distance from the 
primary interaction vertex before scattering, they are 
problematic to associate with the neutrino event. Be¬ 
sides, neutrons typically deposit only part of their energy 
in the detector. Therefore, this assumption may need to 
be revisited in future, when the ongoing experimental 
program |33j [34) brings progress in the understanding of 
detector response to neutrons. 

Very recently, several experimental collaborations have 
started to study in more detail the impact of system¬ 
atic uncertainties on their CP violation sensitivities. In 
Ref. [35], for instance, a fast Monte Carlo was used to 
estimate the impact of the detector performance and nu¬ 
clear effects on the neutrino energy resolution for the 
Deep Underground Neutrino Experiment (DUNE). Nev¬ 
ertheless, when computing the sensitivities to different 
observables, the assumptions which could affect the en¬ 
ergy resolutions were kept fixed in the analysis. Our ap¬ 
proach is different from the one considered in previous 
references: instead, in the present work we consider the 
possible effect in the analysis if the assumptions used 
to get the neutrino energy resolution were very different 
from expectation. 

Imperfect detection capabilities induce a non-vanishing 
probability that an event of a true energy E v ends up 
being reconstructed with a different energy E lec . We en¬ 
code them in a set of migration matrices, calculated as 


3 


in Ref. [31], using genie 2.8.0 [3BJ with the vT modules 
package m- Note that our results inevitably are subject 
to uncertainties coming from nuclear effects. Should the 
argon nucleus be employed as the target, those uncertain¬ 
ties would not be possible to estimate, due to scarcity of 
reported neutrino cross sections. Therefore, in order to 
minimize nuclear uncertainties of our results, we consider 
the carbon target, for which a number of the extracted 
cross sections is available. This allows us to discuss the 
role of detector effects in an unambiguous way. 

The considered experimental setup consists on a wide¬ 
band neutrino beam produced mainly from pion and kaon 
decays, aimed at a 40 kton detector located at a distance 
of L = 1300 km from the source. The neutrino fluxes 
used in this Letter correspond to the 80 GeV beam con¬ 
figuration from Ref. [55], with an assumed beam power 
of 1.08 MW. The background implementation follows 
Ref. [38] as well. No migration matrices are used for the 
background events, which are always smeared according 
to a Gaussian with a(E u ) = 0.15As for the sig¬ 
nal efficiencies, since the detection of neutrino and an¬ 
tineutrino CC events depends on the ability to observe 
and tag only the associated charged lepton, we use the 
same signal efficiencies as in Ref. [35] (80%), where this is 
taken into account. The energy of all particles produced 
in the event (both the charged leptons and the hadrons) 
are then smeared according to a Gaussian, as explained 
in Ref. |31j . before reconstructing the neutrino energy. 
Detection thresholds and efficiencies for all hadrons are 
implemented as well, following Ref. [3T|. The hadron 
thresholds and efficiencies will affect the smearing of the 
events in reconstructed neutrino energy, but not the total 
event rates. 

A total of 6 years (3 in positive horn focus¬ 
ing/neutrino running mode, and 3 in negative horn focus- 
ing/antineutrino running mode) are considered. Under 
these assumptions, the total number of events in the neu¬ 
trino (antineutrino) running modes with reconstructed 
energies between 0.6 GeV and 6 GeV is: 740 (286) signal 
events, 114 (67) intrinsic beam v e and v e events, 67 (33) 
misidentified v tl and events, and 65 (38) neutral cur¬ 
rent events. It should be noted that in the antineutrino 
running mode we consider both v^ —> v e and v fl —>• v e 
as the signal events due to the large contribution to the 
signal from wrong-sign events, which are also sensitive to 
CP violation. In the neutrino running mode, however, 
only the v { , — > v e events are considered as part of the 
signal, since the wrong-sign contribution is negligible. 

A modified version m of globes [591 HO] (General 
Long Baseline Experiment Simulator) is used for the os¬ 
cillation analysis. To determine the confidence regions 
and the significance of the signal, a binned y 2 is con¬ 
structed, following the prescription of Refs. [m mi in]. 
The bin size is set to 100 MeV in reconstructed neutrino 
energy. In this work, however, no near detector is consid¬ 
ered. Instead, we make rather aggressive assumptions for 



FIG. 1: (color online). Reconstructed energy distributions 
obtained for v e deep-inelastic scattering (DIS) events with 
true energy of 2.95 GeV. The distributions neglecting the shift 
due to the missing energy (dot-dashed line), and accounting 
for its 50% (dashed line) are compared to the full calculations 
(solid line). 


the systematic uncertainties, and assume that the near 
detector will be able to achieve these goals. Two sets of 
systematic uncertainties are considered for the signal: a 
normalization (bin-to-bin correlated) and a shape (bin- 
to-bin uncorrelated) uncertainty. A prior at the 2% level 
is considered for both of them, following Ref. [22] 38] . 
As for the background, only a global normalization un¬ 
certainty, at the 5% level, is considered. 

All oscillation parameters are kept fixed in our sensi¬ 
tivity calculations; the conclusions are not expected to 
be qualitatively affected if marginalization over the rest 
of oscillation parameters is performed. Since the atmo¬ 
spheric parameters are fixed to their current best-fit val¬ 
ues, and we are only interested in the 5c p sensitivity, 
there is no need to include and disappearance chan¬ 
nels in our analysis. Therefore, only the results in the v e 
and v e appearance channels are included in our fits. It 
should be kept in mind, though, that the measurement 
of the disappearance parameters may be significantly af¬ 
fected by either an incorrect estimate of nuclear effects 
and/or by an inaccurate detector calibration, as it was 
pointed out in Refs. [sn sb si], among others. An in¬ 
correct determination of the disappearance parameters 
would unavoidably affect the extraction of the value of 
the CP-violating phase from appearance measurements. 

The true event rates are obtained taking into account 
realistic detection capabilities which are implemented us¬ 
ing the migration matrices obtained from Monte Carlo 
events. Therefore, the neutrino energy is not recon¬ 
structed around the true energy but around a lower value 
instead, owing to the energy carried away by unobserved 
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FIG. 2: Effect of an underestimation of the missing energy 
in the calorimetric energy reconstruction on the confidence 
regions in the (#13, 5) plane, see text for details. The true 
values of the oscillation parameters are indicated by the dot, 
and are the same for all contours shown. 


particles in the final state. 

The fitted event rates are smeared using a different 
function. In the ideal case where no particle escapes de¬ 
tection, the neutrino energy would be smeared according 
to a Gaussian distribution centered around the true neu¬ 
trino energy, whose width depends on the energy smear¬ 
ing of the different particles observed. In our analysis, 
the event rates used to fit the data are smeared using 
a linear combination between the two cases described 
above: the realistic scenario where migration matrices 
are used, and the ideal case with a Gaussian smearing 
around the true energy. By varying the coefficients in 
this linear combination, the effective smearing function 
obtained can be deformed smoothly from one situation 
to the other. In this way, we introduce a way to manu¬ 
ally tune the amount of missing energy in the oscillation 
analysis, while at the same time we account for the effect 
of realistic energy resolutions of the detector. 

To illustrate how the energy reconstruction is affected 
by the missing energy, in Fig. |T] we show an example 
for deep-inelastic v e scattering at the true energy E v = 
2.95 GeV. The solid line presents the reconstructed- 
energy distribution calculated from the Monte Carlo sim¬ 
ulations with all detector effects. Should no energy be 
missing, the distribution would be centered at the true 
value of the neutrino energy, as the dot-dashed curve. A 
common way used in the literature to parametrize the 
resolution in neutrino energy in oscillation experiments 


is by using a Gaussian function with a simple function for 
its standard deviation: a{E v ) = a + ^yfEv + r yE v , where 
E„ is the true neutrino energy in GeV. Typical values 
used in phenomenological studies of liquid argon detector 
experiments are <j(E„) = 0.15 y/E^, see e.g. Refs. pCIf 
146j . In our case, we use the migration matrices which 
have been obtained from the event generator, and fit the 
result to a Gaussian with a width in the above form. In 
the case of u e DIS events (i.e., the dot-dashed curve in 
Fig. 0 , the best-fit to the matrices is given by a Gaussian 
with standard deviation a(E v ) = 0.158E U + 0.13 yfE^. 
Finally, the dashed curve, obtained from linear interpo¬ 
lation between the dot-dashed and solid lines, represents 
an intermediate situation in which 50% of the missing en¬ 
ergy is accounted for: the two distributions used in the 
linear interpolation do have the same width, while their 
central value differs due to the impact of missing energy 
in the events. It should also be noted that, for each type 
of neutrino interaction considered in this work, the width 
of the distribution obtained when computing the migra¬ 
tion matrices is generally different. 

Based on Monte Carlo studies, the hadronic energy un¬ 
certainty in the MINOS experiment has been estimated 
not to exceed 8.2% |47| . However, in view of the re¬ 
ported difficulties with the description of nuclear effects 
in modern simulations |48| , our results are presented for 
uncertainties up to 30%. 

The allowed confidence regions from the oscillation 
analysis are shown in the (0 13, S) plane in Fig. [2] In this 
figure, the different contours have been obtained under 
different assumptions regarding the ability of the exper¬ 
iment to determine the missing energy involved in the 
events. The shaded area corresponds to the correct re¬ 
sult, where all the missing energy in the events is per¬ 
fectly estimated in the fit. The solid, dashed, and dot- 
dashed lines represent the results obtained when 90%, 
80%, and 70% of the missing energy is correctly ac¬ 
counted for, respectively. Our results show that even 
a 20% underestimation of the missing energy introduces 
a sizable bias in the extracted Scp value. Should an ex¬ 
perimental analysis suffer from a 30% underestimation 
of the missing energy, it would exclude the true value of 
Scp at a confidence level between 2 and 3er. 

The legend in Fig. [2] also shows the values of the y 2 
for the best fit (#13, S) points divided by the effective 
number of degrees of freedom, i.e., the number of data 
bins minus the number of parameters extracted from the 
data. In an actual experiment, this ratio would give an 
additional contribution to the goodness of fit. A large 
enough contribution would indicate that the model used 
to fit the data is not correct. Our results indicate that 
such contribution would be small enough that, from a 
fit to the far detector data alone, it would be virtually 
impossible to realize that the energy carried away by un¬ 
detected particles is being underestimated in the fit. 

In summary, we have analyzed the impact of missing 
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energy on determination of the CP-violating phase in a 
long-baseline neutrino appearance experiment employing 
the calorimetric method of energy reconstruction. The 
main source of missing energy are neutrons and other 
hadrons escaping detection when realistic detection ca¬ 
pabilities are taken into account. Our results suggest 
that an underestimation of missing energy by as little as 
20% may result in a bias of around 1 standard devia¬ 
tion in the extracted value of Sqp- As a final remark, we 
would like to emphasize that, although the configuration 
considered in our analysis is meant to be similar to the 
design of DUNE, clearly much more detailed studies are 
necessary to draw quantitative conclusions for a specific 
detector setup. In particular, our treatment of the miss¬ 
ing energy uncertainty—assumed to be equal for neutri¬ 
nos and antineutrinos and independent of the energy and 
interaction channel -may be regarded as simplistic. As 
more realistic sensitivity estimates would require an ac¬ 
curate knowledge of the detector response and inclusion 
of nuclear-model uncertainties, out of necessity, we leave 
them for future investigations within experimental col¬ 
laborations. 
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